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Abstract. Computations with a future quantum computer will be imple- 
mented through the operations by elementary quantum gates. It is now well 
known that the collection of 1-bit and 2-bit quantum gates are universal for 
quantum computation, i.e., any n-bit unitary operation can be carried out by 
concatenations of 1-bit and 2-bit elementary quantum gates. 

Three contemporary quantum devices— cavity QED, ion traps and quan- 
tum dots-have been widely regarded as perhaps the most promising candidates 
for the construction of elementary quantum gates. In this paper, we describe 
the physical properties of these devices, and show the mathematical deriva- 
tions based on the interaction of the laser field as control with atoms, ions or 
electron spins, leading to the following: 

(i) the 1-bit unitary rotation gates; and 

(ii) the 2-bit quantum phase gates and the controlled-not gate. 

This paper is aimed at providing a sufficiently self-contained survey account 
of analytical nature for mathematicians, physicists and computer scientists to 
aid interdisciplinary understanding in the research of quantum computation. 



1. Introduction 

The design and construction of the quantum computer is a major project by the 
scientific community of the 21 st Century. This project is interdisciplinary by nature, 
requiring close collaboration between physicists, engineers, computer scientists and 
mathematicians. 

Physicists have taken the lead in the design of quantum devices utilizing the 
peculiar properties of quantum mechanics. Most or perhaps all of such devices will 
eventually become micro- or nano- fabricated. Their controls are achieved through 
the tuning of laser pulses. As hardware, many types of material or devices have 
been either used or proposed: NMR, cavity QED, ion or atom traps, quantum dots, 
SQUID, etc. See some discussions in 20, Chap. 7], |27| . for example. NMR (nuclear 
magnetic resonance) is the first scheme utilized by researchers to demonstrate the 
principle of quantum computing through quantum search. NMR seems to have 
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provided the most successful demonstration of quantum computing so far. However, 
this kind of "quantum computing in a coffee mug" utilizes bulk material, which is 
not considered as really quantum by many people. In addition, NMR has a severe 
weakness in the lack of scalability. Concerning SQUIDs 1 (superconducting quantum 
interference devices), there has been major progress in the design and control of 
such devices recently. However, the authors' analytical knowledge about SQUIDs 
still appears very limited for the time being and, thus, any detailed mathematical 
representations must be deferred to a future account. 

We focus our attention on the following three quantum devices: cavity QED, 
ion traps and quantum dots. They seem to have received the most attention in 
the contemporary literature, and have been widely identified by many as the most 
promising. The writing of this paper all began with our simple desire to try to 
understand, analytically, how and why such devices work. Analytical studies of 
these devices are scattered in many references in the physics literature, which ad- 
mittedly are not easy reading for most mathematicians. Indeed, this is the most 
commonly encountered difficulty in any interdisciplinary research. But even many 
physicists specializing in quantum devices find certain difficulty in their attempt to 
understand the working of other quantum devices, if such devices do not fall exactly 
within their own specialty. Therefore, our objectives in this paper are three-fold: 

(1) Provide a sufficiently self-contained account of the physical description of 
these primary contemporary quantum computing devices and the deriva- 
tion of their mathematical representations. 

(2) Supply the control-theoretic aspect of quantum control via the shaping 
of laser pulses in quantum computing, which is the main theme of the 
conference. 

(3) Write an accessible survey of three important quantum devices for math- 
ematicians and other scientists who are interested in understanding the 
basic interdisciplinary link between physics, mathematics and computer 
science, in a mathematically rigorous way, for quantum computing. 

Our main emphasis will be on the mathematical modeling. As far as the physics 
is concerned in this paper, we acknowledge the many segments that we have ex- 
tracted, directly or indirectly, from the original sources cited in the references. New 
devices and designs emerge almost daily, and we must concede that the list of refer- 
ences is neither exhaustively comprehensive nor totally up-to-date as any attempt 
to achieve such by us is nearly impractical, or even impossible given the time and 
various other constraints. We nevertheless believe that the basic mathematics un- 
derlying most of the quantum computing devices remains largely similar, and hope 
that we have provided enough ideas and information for the readers to trace the 
literature on their own. 



The SQUID can be viewed as a high-Q microwave frequency LC resonator. The first excited 
state, with excitation energy fiw where u) is the microwave frequency, is doubly degenerate in 
zero magnetic field. These two states are the qubit states. To lift the degeneracy, Joscphson 
junction(s) are placed in the current loop. Magnetic field can then be used to tune the levels 
back to degeneracy as required for qubit coupling, etc. The coupling is magnetic, as in mutual 
inductance of two current loops. 

The major breakthrough made last year (2001-02) is the single quantum flux gate for readout. 
It is 100% efficient and does not add decoherence. Expect to see rapid progress in the near future 
as a result of this breakthrough. 

We thank Prof. Philip R. Hemmer for the above communication. 
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In the main body of this paper to follow, there are four sections. Section 2 
provides a quick summary of universality results of quantum gates. Sections 3, 4 
and 5 deal with, respectively, cavity QED, ion traps and quantum dots. In the final 
Section 6, we briefly discuss some issues on laser control in quantum computing as 
the conclusion of this paper. 

2. Universality of Elementary Quantum Gates 

We begin this section by introducing some basic notations used in quantum 
mechanics, as such notations do not appear to be familiar to a majority of mathe- 
maticians. A useful reference can be found in |25j . 

Let TL be a complex Hilbert space with inner product ( , ). For any subset 
S C TL, define 



span S — closure of 




6C,^eS,i = i,2, 



, n, for all ipj G S 



(2.1) 



in TL, 



where in the above, C denotes the set of all complex numbers. Elements in TL , such 
as ipj in (|2.1|l . are called vectors. The inner product of two vectors ipj an d ipk is 
(ipj, ipk)- However, in quantum mechanics, the Dirac bra-ket notation writes the 
vector ipk as |fc), pronounced ket k, where k labels the vector. The role of ipj in the 
inner product (ipj, ipk) is that ipj is identified as a linear functional on TL by 

(2.2) iP): TL^C; ififa) = <^,^fc>- 

The Dirac bra-ket notation writes ip* as (j\, pronounced bra j. The inner product 
is now written as 

(2.3) (iPj,ip k ) = (j\k). 

Let A: TL — > TL be a linear operator. Then in Dirac's notation we write (ipj,Aipk) 
as (j|A|fc). 

The set of all kets (\k)) generates the linear space TL, and the set of all bras 
also generates the dual linear space TL* . These two linear spaces are called, 
respectively, the ket space and the bra space, which are dual to each other. The 
natural isomorphism from H onto H* , or from H* onto H, is called "taking the 
adjoint" by physicists. An expression like (j|^4|/c) has an intended symmetry: one 
can think of A acting on \k) on the right as a matrix times a column vector, or A 
acting on (j| on the left as a row vector times a matrix. 

Any 2-dimensional complex Hilbert space TL is isomorphic to C 2 . The standard 
basis of C 2 is 

(2.4) 

In quantum mechanics, the preferred notation for basis is, respectively, |0) and |1), 
where |0) and |1) generally refer to the "spin" state of a certain quantum system 
under discussion. "Spin up" and "spin down" form a binary quantum alternative 
and, therefore, a quantum bit or, a qubit. These spin states may then further be 
identified with e\ and e2 through, e.g., 





1 




0~ 


ei = 





i e 2 = 


1 



(2.5) 



|0) = 



ID = 
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For a column vector, say ei, we use superscript T to denote its transpose, such as 
ej = [10], for example. The spin states |0) and |1) form a standard orthonormal 
basis for the (single) qubit's ket space. The tensor product of a set of n vectors 
\Zj) G C?, specified by the quantum numbers zj for j = 1,2, ...,n, is written 
interchangeably as 

(2.6) |zi) (g) \z 2 ) ® ••• ® \z n ) = \zi)\z 2 ) ■ ■ ■ \z n ) = \z x z 2 ■ ■■z n ). 

The tensor product space of n copies of C2 is defined to be 



(2.7) 



= span{|jij 2 • • -On) I jk e {0,1}, fc = 1,2,.. . ,n}. 



It is a 2™-dimcnsional complex Hilbert space with the induced inner product from 
C 2 . A quantum state is a vector in (C 2 )®" (or in any Hilbert space Tt) with the 
unit norm 2 . A quantum state is said to be entangled if it is not a tensor product 
of the form l|2.6[l . A quantum operation is a unitary linear transformation on Ti. 
The unitary group U{TL) on TL consists of all unitary linear transformations on TL 
with composition as the (noncommutative) natural multiplication operation of the 
group. 

It is often useful to write the basis vectors \jxj2 • ■ ■ jn) S (C 2 )®", jk & {0, 1}, 
k = 1, 2, . . . , n, as column vectors according to the lexicographic ordering: 



|00 ■ • • 00) 



(2.8) 



-10) 



100 - • -01) 



111... 11) 



|00 - - - 10) 



cf. H2.5fl . Any 2x2 unitary matrix is an admissible quantum operation on one 
qubit. We call it a 1-bit gate. Similarly, we call a 2 k x 2 fe unitary transformation 
a fc-bit gate. Let a 1-bit gate be 



(2.9) 



U 



WOO Uqi 

u w u n \ ' 

we define the operator A m (U) |2] on (to + l)-qubits (with m = 0, 1, 2, . . .) through 

its action on the basis by 

(2.10) 



A m (U){\xix 2 ■ ■ -x m y)) = 



\x\x 2 ■ 



•x m y) 



if A? =1 x k = 0, 



u 0y \xix 2 ■ ■ ■ x m 0) + ui y \xix 2 ■ ■ ■ x m l) if Af =1 x k = 1, 



2 This type of quantum state is actually a more restrictive type, called a pure state. A pure 
state, say ket \k), could be identified with the projection operator \k)(k\. Then a general state is 
a convex sum of such projectors, i.e., Yl c k\k)(k\ with c k > and c fe = 1. 
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where "A" denotes the Boolean operator AND. The matrix representation for 
A m (U) according to the ordered basis l|2.8|) is 

1 



(2.11) 



A m (C/) = 



o 



1 



o 



uoo 



"01 



2'™+ 1 x2 m + 1 



This unitary operator is called the m-bit-controllcd U operation. 
An important logic operation on one qubit is the NOT-gate 



(2.12) 



cf. the Pauli matrices in (|5.2|l of Section From a x , define the important 2- 
bit operation Ai (cr x ), which is the controlled-not gate, henceforth acronymed the 
CNOT gate. 

For any 1-bit gate A, denote by A(j) the operation (defined through its action 
on the basis) on the slot j: 
(2.13) 



A{j)\xix 2 ■ 



*2 



<£> \Xj 



[A\ Xj 



<8> \x r , 



Similarly, for a 2-bit gate B, we define B(j, k) to be the operation on the two qubit 
slots j and k |4j . A simple 2-bit gate is the swapping gate 



(2.14) 



U sw \xix 2 ) = \x2X1), for xi,x 2 G {0, 1}. 



Then U svr (j, k) swaps the qubits between the j-th and the fc-th slots. 

The tensor product of two 1-bit quantum gates S and T is defined through 

(2.15) {S®T)(\x)®\y) = (S\x))®(T\y)), for x, y G {0, 1}. 
A 2-bit quantum gate V is said to be primitive ([4] Theorem 4.2]) if 

(2.16) V = S(g>T or V = (S ® T)U SW 

for some 1-bit gates S and T. Otherwise, V is said to be imprimitive. 

It is possible to factor an n-bit quantum gate U as the composition of fc-bit 
quantum gates for k < n. The earliest result of this universality study was given 
by Deutsch 8 in 1989 for k = 3. Then in 1995, Barenco 2 , Deutsch, Barenco 
and Ekert |2J, DiVincenzo |12| and Lloyd [16] gave the result for k = 2. Here, we 
quote the elegant mathematical treatment and results by J.-L. Brylinski and R. K. 
Brylinski in 0]. 

Definition 2.1 (0 Def. 4.2]). A collection of 1-bit gates A, and 2-bit gates 
Bj is called (exactly) universal if, for each n > 2, any n-bit gate can be obtained 
by composition of the gates Ai(£) and Bj(£, m), for 1 < £, m < n. □ 

Theorem 2.1 (g| Theorem 4.1]). Let V be a given 2-bit gate. Then the fol- 
lowing are equivalent: 

(i) The collection of all 1-bit gates A together with V is universal. 

(ii) V is imprimitive. □ 
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A common type of 1-bit gate from AMO (atomic, molecular and optical) devices 
is the unitary rotation gate 



-ie^ sin( 



-te '^sm( 
cos 9 



< 0, (ft < 2tt. 



(2.17) Ue,* = 

We have (the determinant) det Ug^ = 1 for any 9 and (ft an( h thus, the collection 
of all such t/#,0 is not dense in U(C 2 ). Any unitary matrix with determinant equal 
to 1 is said to be special unitary and the collection of all special unitary matrices 



on 



in , SU({C 2 )® n ), is a proper subgroup of U((C 



2\®ra\ 



It is known that the 



gates Z7e,0 generate SU(C 2 ). 

Another common type of 2-bit gate we shall encounter in the next two sections 
is the quantum phase gate (QPG) 



(2.18) 



1 

1 

1 









e «7 



< t) < 2tt. 



Theorem 2.2 Theorem 4.4]). The collection of all the 1-bit gates Uo,<j>, 
< 0, (ft: 4> 5~ 27r, together with any 2-bit gate Q v where -q ^ (mod 2tt), is 
universal. □ 

Note that the CNOT gate Ai(a x ) can be written as 

(2.19) Ai(o- x ) = C/ 7r/ 4, V 2(2)Q.C/, / 4,-./ 2 (2); 

cf. [101 (2.2)]. Thus, we also have the following. 

Corollary 2.3. The collection of all the 1-bit gates Ug^, < 9, (ft < 2n, 
together with the CNOT gate Ai(cr x ), is universal. □ 

Diagramatically, the quantum circuit for (|2.19|) is given in Fig. 12.11 





Qit 


















Un it 

4 ' 2 






Uit It 

4 ' 2 













FIGURE 2.1. The quantum circuits for equation l|2.19|) . The top 
wire represents the leading qubit (which is the most significant), 
while the second wire represents the second qubit. 



3. Two-Level Atoms and Cavity QED 

The main objective of this section is to show that the 1-bit unitary gates (|2.17() 
and the 2-bit QPG l|2.18|l can be implemented using the methods of, respectively, 
2-level atoms and cavity QED. For clarity, we divide the discussions in three sub- 
sections. 
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3.1. Two-Level Atoms. The atomic energy levels are very susceptible to 
excitation by electromagnetic radiation. The structure of electronic eigenstates 
and the interaction between electrons and photons can be quite complicated. But 
under certain assumptions we have in effect a 2-level atom. These assumptions are 
as follows: 

(i) the difference in energy levels of the electrons matches the energy of the 
incident photon; 

(ii) the symmetries of the atomic structure and the resulting "selection rules" 
allow the transition of the electrons between the two levels; and 

(iii) all the other levels are sufficiently "detuned" in frequency separation with 
respect to the frequency of the incident field such that there is no transition 
to those levels. 

Then the model of a 2-level atom provides a good approximation. 

II) 



transition 
frequency 



10) 

Figure 3.1. A 2-level atom, where a bound electron in an atom 
interacts with an external electromagnetic field. 

The wave function ip(r, t) of the electron of the 2-level atom, which is coupled 
to the external electric field E(r, t) by its electric charge e, as shown in Fig. I3.ll 
obeys the Schrodinger equation 

(3.1) ih §i^ *) = H ^ r > ^ = {Ho + Hl M r > *)' 

where \ip(r, t)\ 2 is the probability density of finding that electron at position r = 
(x, y, z) at time t, and 

(3.2) Hq = V 2 + V(r) is the differential operator for the unperturbed 

^ TO Hamiltonian of the electron, with the gradient 

operator V = (d/dx, d/dy, d/dz) and the elec- 
trostatic potential V(r) between the electron 
and the nucleus, 

(3.3) Hi = — er ■ E(ro, t) is the interaction Hamiltonian between the 

field and the electron of the atom. 
Remark 3.1. The wavelength of visible light, typical for atomic transitions, 
is about a few thousand times the diameter of an atom. Therefore, there is no 
significant spatial variation of the electric field across an atom and so we can re- 
place E(r,t) by E(ro,t), the field at a reference point inside the atom such as the 
position of the nucleus or the center of mass. Consistent with this long- wavelength 
approximation (known as dipole approximation in the physics literature) is that the 
magnetic field B satisfies B(r,t) = 0, so that the Lorentz force of the radiation field 



£cos (Dt)c x 



monochromatic 
electromagnetic field 
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on the electron is eE(ro,t), the negative gradient of — er ■ E(ro,t), which adds to 
the potential energy V(r) in the Hamiltonian operator. □ 

Let the external electromagnetic field be a monochromatic plane-wave field 
linearly polarized along the x-axis, which interacts with the atom placed at ro = 0. 
The electric field takes the form 

(3.4) E(0, t)=£ cos(vt)e x ; e x = [1,0, 0] T . 

Let |1) and |0) represent upper and lower level states of the atom, which are eigen- 
states of the unperturbed part Hq of the Hamiltonian corresponding to the eigen- 
values huji and fajJo, respectively, with lo = u>i — u)q, cf. Fig. 13.11 

For the unperturbed 2-level atom system, i.e., without the external electro- 
magnetic field, let the eigenstates of the electron be \tpo) — |0) and \ipi) = |1), i.e., 



(3.5) 



Hotyo) = E \fa), or H \0) = fiw o |0); <0|0) = 1, 
H \fa)=E 1 \fa), OT Wi) = ^i|l); <1]1> = 1. 



We invoke the assumption of 2-level atom that the electron lives only on states |0) 
and |1): 

(3.6) |^(*)> = C (*)|0> + C 1 (t)|l), 

where Co(t) and C\[t) are complex- valued functions such that |Co (t) | 2 + |Ci(i)| 2 = 
1, with | C (i) | 2 and |Ci(£)| 2 being the probability of the electron in, respectively, 
state |0) and |1) at time t. 

Since Hi (after quantization) is an operator effective only on the subspace 

(3.7) S = span{|0),|l)}, 

Hi is the zero operator on S 1 - (the orthogonal complement of S in the Hilbert space 
7i which has an orthonormal basis consisting of all eigenstates of Hq). We utilize 
the property that on S the projection operator is |0)(0| + |1)(1|, and obtain 

H\ = — er ■ E(0, t) = —e(xe x + ye y + ze z ) ■ {£ cos(vt)e x ) 
= -e£xcos(vt) (by (221 and J33J0 
= - e f[|0)(Q| + |l){lH|0)(0| + |l)(l|]co S (W) 
= -e£{[(0N0)]|0)(0| + [(0Nl)]|0)(l| 

(3.8) +[(l|a:|0)]|l)(0| + [(lHl)]|l)(l|}co S (^). 

But from the symmetry property and the selection rules that ^(f)! 2 = \4>j(~ r )\ 2 
for j = and 1, we have 



oo oo oo 



(3.9) (0|x|0) = / / / x\ip (r)\ 2 dxdydz = 



— oo — oo — oo 



since x = r-e x is an odd function of r whereas li/'ol'")! 2 is even. Similarly, (l|x|l) = 0. 
Now, write Poi = e (0|x|l) and P w =e(l|a;|0). Then P w = P ox , To solve the partial 
differential differential equation 1)3. 1|) . we need only focus our attention on the 2- 
dimensional invariant subspace S and, thus, a dramatic reduction of dimensionality. 
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Substituting Q3.6|l into (|3.1|) . utilizing l|3.5|) and Ij3.8(l and simplifying, we obtain a 
2x2 first order linear ODE 



(3.10) 

where 

(3.11) 

and 

(3.12) 



d 


~C (t) 




di 


C 1 {t)_ 









-tu> 
iH,Re 1 ^ cos vt 



i£lne *^ cos vt 
—iu>i 



C (t) 



\Pio\£ 



= the Rabi frequency, 



phase of P w : P 1Q = \P w \e 



Equation i|3.10fl has time- varying coefficients and exact solutions of such equations 
are not always easy to come by. Here, we make another round of approximation, 
called the rotating wave approximation, by first setting 



(3.13) 



co(t) = C (t)e 



^o* ci(t) = C\{t)e luJlt , 



and substituting them into (|3.1Q(I and then dropping terms involving e ± l ("+' y ) f ) 
where ui = ujq — w\ is called the atomic transition frequency. (From an experimen- 
tal point of view, such e ±l ( UJ + u ) t terms represent high frequency oscillations which 
cannot be observed in laboratory instrumentation. Furthermore, integrals involving 
highly oscillatory terms make insignificant contributions. In view of the fact that 
other non-resonant terms have already been discarded, these e ±l ( ul + u ) t terms rep- 
resenting highly non-resonant contributions should be neglected for consistency.) 
Then we obtain 



(3.14) 



ci(*) = -^-< 



i4> p —i{ui—v)t 



co(t), 



which in turn leads to a second order single ODE with constant coefficients 
(3.15) 



O 2 
"it 



c - i(u - v)cq + —^co = 0. 

Solving for cq and c\ in l|3.14|l and (|3.15|) . we obtain the explicit solution in terms 
of the initial condition (co(0), ci(0)): 
(3.16) 



co(t) 



where Q = 
(3.17) 
and, thus 



e 1 ^* [cos (U 



ft . j gi< 



Q 

' sin 



in (ft)] 



4 



in (ft) 
[cos (§*)+*£ sin (§t)] 



co(0) 
ci(0) 



y/ttjt vf 



At resonance, 

A = uj - v = 0, 



(3.18) 

The above yields 
(3.19) 



A 
9. 



0, Oj 



n. 



c (4)" 




cos (ft) 


_ci(4)_ 




ie^sin(f 



te 



sin (§4) 



cos (§4)" 





'co(O)' 




ci(0)_ 
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(Note that if, instead, we relate Cj(t) to Cj(0) for j = 0, 1, then because of lj3.13|l . 
unlike the (special) unitary matrix in (|3.19|1 whose determine is 1, the corresponding 
unitary matrix will contain some additional phase factor (s).) Write 



(3.20) 



and re-name 6' as 



Then the above matrix becomes 



(3.21) 



cos 9 



-ie v smt 
cos 9 



which is the 1-bit rotation unitary gate. 
Note that 9 in l|3.21|l depends on t. 



Remark 3.2. If we write the Hamiltonian as ([210 (5.2.44), p. 157]) 

Ml, 



(3.22) 



then we obtain 



(3.23) 



H 



-**|1><0| + e**|0><l| 



cos 
ie 1 '*' sin 



an(Ut) 



ie 



sin 



cos (■ 



which is the same as the matrix in l|3.19|l . This Hamiltonian in (|3.22l) is thus 
called the effective (interaction) Hamiltonian for the 2-level atom. It represents the 
essence of the original Hamiltonian in (|3.1|l after simplifying assumptions. From 
now on for the many physical systems under discussion, we will simply use the 
effective Hamiltonians for such systems as supported by theory and/or experiments, 
rather than to derive the Hamiltonians from scratch based on the Schrodinger's 
equations coupled with the electromagnetic field. □ 

3.2. Quantization of the Electromagnetic Field. The quantization of the 
electromagnetic radiation field as simple harmonic oscillators is important in quan- 
tum optics. This fundamental contribution is due to Dirac. Here we provide a 
motivation by following the approach of [251 Chap. 1, pp. 3-4]. We begin with 
the classical description of the field based on Maxwell's equations. These equations 
relate the electric and magnetic field vectors E and B, respectively. Maxwell's 
equations lead to the following wave equation for the electric field: 



(3.24) 



1 d 2 E 



0. 



along with a corresponding wave equation for the magnetic field. The electric field 
has the spatial dependence appropriate for a cavity resonator of length L. We 
take the electric field to be linearly polarized in the a;-direction and expand in the 
normal modes (so-called in the sense that they constitute orthogonal coordinates 
for oscillations) of the cavity 



(3.25) 



E x (z, t) = \] AjQj{t) sm(kjz), 
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where qj is the normal mode amplitude with the dimension of a length, kj = jir/L, 
with j = 1,2,3,.. ., and 



/2^ 2 mA 1/2 

(3.26) A ^[-^r) ' 

with Vj = jwc/L being the cavity eigenfrequency, V = LA (A is the transverse area 
of the optical resonator) is the volume of the resonator and rrij is a constant with 
the dimension of mass. The constant rrij has been included only to establish the 
analogy between the dynamical problem of a single mode of the electromagnetic field 
and that of the simple harmonic oscillator. The equivalent mechanical oscillator 
will have a mass rrij , and a Cartesian coordinate qj . The nonvanishing component 
of the magnetic field B y in the cavity is obtained from Eq. I|3.25|l : 

(3.27) B^X^feWM)- 



The classical Hamiltonian for the field is 

(3.28) H * = \l dT ^ E l + — B l 

2 Jv 



where the integration is over the volume of the cavity. It follows, on substituting 
from (|3.25|) and (|3.27[) for E x and B y , respectively, in (|3.28|l . that 



(3.29) 



2 ( m ^J1] 



2„2 



where pj = rrijqj is the canonical momentum of the jth mode. Equation (|3.29|) 
expresses the Hamiltonian of the radiation field as a sum of independent oscilla- 
tor energies. This suggests correctly that each mode of the field is dynamically 
equivalent to a mechanical harmonic oscillator. 

The present dynamical problem can be quantized by identifying qj and pj as 
operators which obey the commutation relations 

(3.30) ki,qj'] = \Pj,Pj'}=0- 

In the following we shall restrict ourselves to a single mode of the radiation 
field modeled by a simple harmonic oscillator. The corresponding Hamiltonian is 
therefore given by 



(3.31) 



it 1 2 i 1 2 2 

2m 2 
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where 



h d 

p = the particle momentum operator = — — , 

i dx 



m = a constant with the dimensions of mass, 

x — the position operator (corresponding to the q variable above), 
v = the natural (circular) frequency of the oscillator 
(3.32) (a parameter related to the potential depth). 

So the eigenstates ip of the Schrodinger equation satisfy 

h 2 d 2 ip 1 



(3.33) 



mv 2 x 2 il) = Eip. 



2m dx 2 

Let us make a change of variables 
(3.34) y = 

then H3.33fl becomes 
(3.35) 

We now define two operators 
(3.36) 



E 



^ d 2 ip 
dy 2 



y ip 



-\ip. 



a= -j=(d/dy + y), a) 



V2 



(d/dy-y). 



Note that cr is the Hermitian adjoint operator of a with respect to the L 2 (R) inner 
product. Then it is easy to check that for any sufficiently smooth function <f> on K, 



(3.37) 



(oft* 



a) a) 6 = 



(3.38) 



the commutator of a with a) — [a, a'] = 1, 



where 1 denotes the identity operator but is often simply written as (the scalar) 1 
in the physics literature. It is now possible to verify the following: 
(i) Let£=-i(^-y 2 ). Then 

(3.39) H = <Ja+\l = aa) 



(ii) If ip\ is an eigenstate of H satisfying 

(3.40) Htj)x = A?/>a (i-e., 

then so are ai[>\ and <vip\: 

H(ai> x ) = (\-l)aip x , 
H{aH x ) = (X + l)a^ k . 

(iii) If Aq is the lowest eigenvalue (or energy level) of H, then 



(3.41) 
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(iv) All of the eigenvalues of H are given by 




(v) Denote the n-th eigenstate ip n by \n), for n = 0, 1, 2, . . . . Then 



(3.43) |n) = 7^10). 




Furthermore, 

a^a|n) = n\n), n = 0, 1,2, . . . ; 



t| n ) — y/n + 1 \n+l), 7i = 0,l,2,...; 
(3.44) a\n) =Vn\n-l), n = l,2,...; a|0) = 0. 

(vi) The completeness relation is 



oc 



(3.45) 



l = £|n)R 



n=0 



(vii) The wave functions are 



(3.46) 



1>j(y) = N j H j (y)e 



y = 




j = 0, 1, 2, . . ., where Hj(y) are the Hermite polynomials of degree j, and Nj is a 
normalization factor, 



back to the cc-coordinate from the y-coordinate in l|3.42Jl ) may be inter- 
preted as the presence of n quanta or photons of energy hv. The eigen- 
states | n) are called the Fock states or the photon number states, 
(ii) The energy level E = \fw (from l|3.34|l and (|3.41|) ) is called the ground 
state energy. 

(hi) Because of (I3.44|) . we call a and a', respectively, the annihilation and 
creation operators. □ 

3.3. Cavity QED for the Quantum Phase Gate. Cavity QED (quantum 
electrodynamics) is a system that enables the coupling of single atoms to only a 
few photons in a resonant cavity. It is realized by applying a large laser electric 
field in a narrow band of frequencies within a small Fabry-Perot cavity consisting 
of highly reflective mirrors. See a primitive drawing in Fig. 13.21 

A 3-level atom is injected into the cavity. An electron in the atom has three 
levels, \a), \(3) and I7), as shown in Fig. 13.31 Actually, the state \a) will be used 
only as an auxiliary level because later we will define the states |/3) and I7) as 
the first qubit, |1) and |0), respectively. Once the atom enters the cavity, the 
strong electromagnetic field of the privileged cavity mode causes transitions of the 
electron between \a) and \0), and a photon or photons are released or absorbed in 
this process. Of the photon states |0), |1), |2), . . . inside the cavity only |0) and |1) 
will be important (which physically mean or 1 photo inside the cavity), and they 
define the second qubit. 

The creation and annihilation operators a) and a act on the photon states |n) 
for n — 0, 1, 2, ... , 00, according to l|3.44|l . 



Remark 3.3. 




converting 
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-highly reflective mirror 



Figure 3.2. A Fabry-Perot cavity for electromagnetic radiation. 
Inside the cavity, the mode structure of the radiation field is dras- 
tically changed. In particular, there is a single privileged cavity 
mode that is resonant with the atomic transition and dominates 
the atom-field dynamics so completely that the influence of all 
other modes is negligible. 



V =(0y a + A 



CO 



ap 



IY> 
|a> 

IP> 



cavity 



atom 



Figure 3.3. Diagram for cavity QED. 



The Hamiltonian for the atom-cavity field interaction is given by 

(3.47) H = H Q + H 1 +H 2 , 
where 

Hq = — ^-{\a){a\ — |/3)(/3|) = the atom's Hamiltonian, 

(3.48) Hi = hva) a = the Hamiltonian of the laser electric field of the cavity, 
and 

(3.49) H-2 = Hg(\a)(/3\a + \f3)(a\<v) = the interaction Hamiltonian of 

the laser field with the atom; with g > 0. 

See the derivation of l|3.47|l in |21| . Note that the operators a and a' appearing in 
(I3.47|l operate only on the second bit, while the rest of the operators operate only on 
the first bit. Also, operator \a}(a\ — \f3}(/3\ can be written as a z ; see equation (|5.2|) 
later, because its matrix representation with respect to the ordered basis {\a}, 
is exactly a z . 
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Lemma 3.1. Let the underlying Hilbert space be 

(3.50) H = span{\a, n), |/3, n), \ j, n) \n = 0, 1, 2, . . .} 

Then the Hamiltonian operator i|3.47fl has a family of 2- dimensional invariant sub- 
spaces 

(3.51) V n = span{\a, n — 1), |/3, n)}, n=l,2, ... . 

Indeed, with respect to the ordered basis in l|3.51fl . the matrix representation of H 
on V n is given by 

\ujap + v(n - 1) g^/n 



(3.52) 



H 



jW a /3 + vn 



n = 1,2,3, 



Proof. Even though the verification is straightforward, let us provide some 
details for the ease of future reference. We have, from l|3.44ll . and 13.47fl - H3.49fl . 

H\a,n - 1) = H \a,n - 1) + Hi\a,n - 1) + H 2 \a,n - 1) 
hw a p 



2 

(3.53) = A 
Similarly, 

(3.54) H\f3,n) = h 
Therefore, we obtain (|3.52fl 



\a, n — 1) 



[hv{n — \)\a, n — 1)] + [hg^/n\/3, n}] 



+ v(n— 1) \a,n— 1) + h[g-^n\\(3,n). 



|/3, n) + h{gy/n)\a, n — 1). 



□ 



Lemma 3.2. On the invariant subspace V n in (|3.51fl . the Hamiltonian operator 
H has two eigenstates 



(3.55) 

where 



cos9 n \a,n — 1) — sm9 n \f3,n), 
sm9 n \a, n — 1) + cos 



sin 0„ 



- 2 1. 



COS 



2g^/n 



D 1 " D ' 

D= [(r»„-A) 2 +4 5 2 n] 1 / 2 , 

fi„^(A 2 +4.g 2 n) 1 / 2 , 
(3.56) A = ^ — u; a/ g = the detuning frequency, 

with eigenvalues (i.e., energy levels) 
(3.57) 
such that 

(3.58) H\+) n = £+(„)]+>„, #|->„ = S_(„)|-)„. 

Proof. These eigvenvalues and eigenvectors can be computed in a straight- 
forward way, using the 2x2 matrix representation of H on V n , 13.52fi . with respect 
to the ordered basis of V n chosen as in (I3.51fl . □ 



E±( n ) — h~ 



nv + \^{-v =F fi„) 
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Remark 3.4. The states |+)„ and |— )„ in 1|3.55|1 are called the dressed states 
in the sense that atoms are dressed by electromagnetic fields. These two states are 
related to the splitting of the spectral lines due to the electric field. □ 

Now, let us assume large detuning: 
(3.59) |A|>2 5% AI. 

Also, |A| remains small in comparison to the transition frequency so that approx- 
imate resonance is maintained. Accordingly, this is satisfied only for small values 
of ri. Then 

^ = (A 2+ Vn) 1/2 = Afl + ^V /2 -A + ^; 



(3.60) sinl? n = (fJ„-A)/.D«0; cos0 n = 2g^/n/D m 1. 
Thus, from 

(3.61) \+) n n\a,n-l), |->„ « \/3, n). 



Remark 3.5. By the assumption i|3.58[) of large detuning, from i|3.57|l using 
(|3.56|l and Ij3.59|l . we now have 



E. 



+(«) 



(3.62) 

Similarly, we have 

(3.63) 
such that 

(3.64) 



; a 



n — — \ v 



2g 2 n 



g 2 n 



n - r f - ~{v-uj al3 ) - — 



-^■ + v(n-l) 



E_ 



in) 



vn 



hg 2 n 
A ' 

hg 2 n 



+ 



H\ 



H\a,n- 1) = E. 



+(n) 



A ' 
a, n — 1) 



H\-) n ^H\/3,n)=E_ (n) \0,n), 

where E±(n) are now given by (|3.62|1 and (|3.63(1 . 

Thus, we see that the assumption H3.58JI of large detuning causes the off- 
diagonal terms in the matrix H3.52|l (which are contributed by the interaction 
Hamiltonian H2 in J23EJ0 to disappear. The de facto or, effective interaction 
Hamiltonian, has now become 

hg 2 



(3.65) 



Ho 



A 



•(aa t |a)(a| — a. t a|/3)</3| ) , 



such that H = Hq + H\ + H2 now admits a diagonal matrix representation 



(3.66) 



H\ 



E 



+(«) 







E_ 



(«). 
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with respect to the ordered basis {|a, n — 1), |/3,n)} of V n , with £?±( n ) here given 
by and dHJ. □ 

We now define the first qubit by 

(3.67) |/3) = |1), |7) = |0>, 

where I7) is the state in Fig. 13.31 that is detached (or off-resonance) from the 
Hamiltonian H, i.e., 

(3.68) H\<y) = 0. 

For the second qubit, we choose n = 1 for V n in l|3.51|l and then see that the 
second bit is only or 1, i.e., the cavity can have or 1 photon. The cavity field is 
nearly resonant with the \a) <-> \j3) transition: v = oj a p + A. 

Theorem 3.3. Let the effective total Hamiltonian be 

H = H + H 1 +H 2 = h 



^(| a )<a|-|WI) Wa 



(3.69) -'-^-(aa)\a){a\ - o)a\f3){(3\ 
Then 

(3.70) V- = «po»*{|0, 0>, |0, 1>, |1, 0>, (1, 1» 

is an invariant subspace of H in the underlying Hilbert space <|3.50[1 , such that 

iir|o,o) = 0, #|0,1) = 0, #11,0) = 



(3.71) 



£r|l,l) = £_m|l,l), where E_ (1) = H[-^ + v + ^- 



Consequently, with respect to the ordered basis ofV in <|3.70[l . H admits a diagonal 
matrix representation 



(3.72) H\ 
with the evolution operator 









O 



O 



E -w. 



(3.73) 



iHt/H\ 



o 



o 



exp I — i 



Proof. We have 

(3.74) #|0,0) = H\-y,0) = and H\0, 1) = H\-y, 1) = 0, by 
Also, 

(3.75) g|l,0) = H\3.0) = 0. bv (l3~H6l . 

where V n is chosen to be Vb by setting n = therein. But 

(3.76) H\l,l)=H\0,l) = E_, 1) \l,l) 
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by H3.66(l where V n is chosen to be V\ by setting n 

inzznii follow. 



1 therein. So both (pTf^l and 

□ 



The unitary operator in H3.73|) gives us the quantum phase gate 

"1 



(3.77) 

with 
(3.78) 



1 

o 



1 



o 

P i<t> 



We can now invoke Theorem 12 .21 in Section [5] to conclude the following. 

Theorem 3.4. The collection of 1-bit gates Ug^ in Ij3.21|l from the 2-level 
atoms and the 2-bit gates Q v in H3.77(l from cavity QED is universal for quantum 
computation. □ 



4. Ion Traps 

Ion traps utilize charged atoms with internal electronic states to represent 
qubits. Ion traps isolate and confine small numbers of charged atoms by means 
of electromagnetic fields. The ions are then cooled using laser beams until their 
kinetic energy is much smaller than the inter-ionic potential energy. Under these 
conditions, the ions form a regular array in the trap. Laser beams can be tuned to 
excite the electronic states of particular ions and to couple these internal states to 
the center-of-mass (CM) vibrational motion of the ion array. This coupling provides 
entanglement for quantum computation. 

A Paul, or radio-frequency, ion trap confines ions by the time average of an 
oscillating force F(r, t), where r = (x, y, z), arising from a non-uniform electric field 
E(r) cos(fit) acting on the ion charge e with mass m. In one dimension, (F(X, t)) — 
—d(U(X,t))/dX, where (U(X,t)) is an effective time-averaged potential energy, 
and U(X,t) = e 2 E 2 (X) cos 2 (fit)/(mf2 2 ). The time average is denoted by ( ) and 

here, it is defined by (y(t)) — ^ f Q n y(t)dt, where y(t) is periodic with period 
2-k/Q,. X describes the (relatively slow) motion of a guiding center, about which 
the ions execute small oscillations at frequency f2. The trap is constructed so 
that the guiding center motion in two dimensions (say x and y) is harmonic, with 
frequencies v x and v y which are much smaller than the frequency Q. The trap 
electrodes may confine the ions in a circle or other closed 2-dimensional geometry, 
or a linear trap may be produced by using a dc potential to confine the ions in the 
third spatial dimension z; see Fig. 14.11 

Once confined in the trap, the ions may be cooled using laser light, so that 
they condense into a linear array along the trap axis. Each ion oscillates with small 
amplitude about the zero of the time-averaged radial potential. The electrostatic 
repulsion of the ions along the linear axis of the trap combats a spatially varying 
dc potential applied for axial confinement, so that below milliKelvin (1 x 10~ 3 K) 
temperatures, the ions in the array are separated by several micrometers (lCU 6 ™). 
Under these conditions, the ion confinement time is very long, many hours or days, 
and each ion can be separately excited by a focused laser beam, so that long-lived 
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Figure 4.1. Diagram of a linear ion trap, showing the radio- 
frequency voltage connections on the four rods of the trap, and 
the dc voltage connection for an axial dc potential well. One rod is 
partially cut back to reveal 6 cold ions (hiding behind, not inside, 
the rod) confined in a linear array within the trap. The ions would 
be visible due to scattered light, if they were all simultaneously ex- 
cited by a laser beam tuned to an internal resonance transition, at 
a wavelength in the visible region of the spectrum. The separation 
of the trap rods is on the order of a millimeter. 

internal levels of each ion can serve as the two states of a qubit (Fig. 14. 2|) . Re- 
garding the axial kinetic motion, the ions oscillate axially as a group, so collective 
oscillations such as the CM oscillations couple the individual ions, or qubits, per- 
mitting the implementation of multi-qubit gates based on quantum entanglement 
of the CM motion with each laser-excited qubit state. There are many CM modes 
of oscillation, depending on the oscillations of individual ions or groups of ions, 
each with a characteristic frequency wcm- 

A prototype scheme was proposed by Cirac and Zoller |7|, which is based on 
a confined ion system as outlined above. Assume the quantized axial CM normal 
mode oscillation has frequency wqm and discrete levels nhu>cM with quantum num- 
ber n = 0, 1, 2, . . . (the ground state zero-point energy ?iwcm/2 in Remark l3.3f ii) 
is disregarded) ; see Fig. I4.3f bl . A quantum of oscillation Hujcm is called a phonon. 
Each ion is individually excited by focused standing wave laser beams. Two inter- 
nal hyperfme structure states of the ion (i.e., states coupling the electrons and the 
nucleus of the ion) in the ground term (or lowest electronic set of states), separated 
by the frequency luq, can be identified as the qubit states |0) and |1). Alternatively, 
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|oc) 



|cc) is the auxiliary level 



(D 



10) 



Figure 4.2. Diagram of the energy levels of an idealized ion, 
which might serve as the two qubit levels, plus an auxiliary level 
used in a particular implementation of a quantum circuit. The ac- 
tual levels of a real ion are much more complex, but approximate 
the ideal case when the exciting lasers are properly tuned. 



a ground level and a metastable excited level might be used. By thus coupling only 
these selected states of the ion by electromagnetic fields we are led to a 2-level atom 
analogously as in Section These states can be coupled using electric dipole tran- 
sitions to higher-lying states, and there is also an internal auxiliary level, denoted 
| a), which can be coupled to the qubit states; see Fig. 14.21 Consequently a linear 
iV-ion system can be specified in terms of quantum numbers by the wavefunction 
\ip)N = \h)\h) ■ ■ ■ \]N)\k)cM, where states (of each ion) jx, j 2 , . .. Jn € {0, 1} and 
k = 0, 1, 2, . . . , N — 1 describes the populated quantum state of the CM motion. 
The system is initialized with k = when the ions are very cold, and with all 
qubit states Ji , J2> • • • j Jn set to by inducing appropriate laser transitions in the 
individual ions. 

Fig. 4.3(b) gives a diagram of the harmonic potential well of the axial ion center- 
of-mass vibration, showing the lowest energy levels. The six ions in Fig. 14. 3f a~) 
would populate one of these energy levels, acting as a single entity with the CM 
coordinates. 

The acting laser beam tuned to frequency too causes a dipole transition coupling 
the two qubit levels |0) and |1) in Fig. 14.21 Using the same rotating wave approxi- 
mation as in Sectional we obtain the 1-bit rotation matrix Ug^ as in 1(3.21(1 . 

Therefore, according to the theory of universal quantum computing (Theorem 
12.1(1 in Sectional now all we need to show is an entanglement operation between 
two qubits. 

We assume that two ions corresponding to the two qubits are coupled in the 
Lamb-Dicke regime, i.e., the amplitude of the ion oscillation in the ion trap potential 
is small compared to the wavelength of the incident laser beam. This is expressed 
in terms of the Lamb-Dicke criterion: 



(4.1) 



r]{n/2ujcM) < i 
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Figure 4.3. (a) A diagram showing 6 cold ions in a linear array, 
the arrows show two of the ions individually excited by focused 
laser beam standing waves. One laser is tuned to the internal 
qubit transition frequency loq, while the other laser is tuned to the 
first red sideband frequency loq — wcm , which couples the internal 
transition to the center-of-mass vibration of the ion array at wcm 
in the ion trap. 

(b) A diagram of the harmonic potential well of the axial ion center- 
of-mass vibration, showing the lowest energy levels. The six ions 
in (a) would populate one of these energy levels, acting as a single 
entity with the CM coordinates. 



where 



the Lamb-Dicke parameter, 
a "zero-point" oscillation amplitude, 
at the lowest energy of CM vibration, 
the wave number of the laser radiation, 
the frequency of the CM motion referred to earlier, 
the Rabi frequency (cf. (|3.11l) '). 

When (|4.1|l is satisfied, the "sideband" frequencies u>o ± u>cm arising from the 
CM oscillation in the laser standing wave are resolved (see Fig. 4.4 and the captions 
therein), and can be separately excited by the laser beam; see Fig. I4.3f a'l. Then 
the interaction Hamiltonian coupling the internal qubit states of each ion j, j = 



(4.2) 



= kz 
zo 

2tt/A 

WCM 

n 
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Figure 4.4. oj — wcm is the red sideband frequency which may 
be excited to entangle the qubit with the CM motion. The right 
figure indicates that when the Lamb-Dicke criterion is satisfied the 
sideband frequencies luo ± wcm and u>q, which have finite widths 
due to the Heisenberg uncertainty principle, are resolved, enabling 
separate excitations by the laser beam. (The frequency luq + wcm 
is called the blue sideband frequency; see [28].) 



1, . . . , N, to the CM motion is similar to 1(3.49(1 . with coupling constant g = 2 ^L- : 

(4.3) H 3 = ^Ol^e-^IO^liyM, j = l,..-,N, 

l\ iv 

where a and are, respectively, the annihilation and creation operators (see Sub- 
section [!OJ for the phonons defined by 

a|n)cM = \fn \n - 1)cm, for n = 1, 2, . . . ; a|0)cM = 0|0)cm = 0, 

(4.4) a f |n)cM = Vn+l \n+ 1)cm, for n = 0, 1, 2, . . . ; cf. EH) , 
and 

4> = the angle in the spherical coordinate system (r, 9, <fi) with respect 

(4.5) to the rr-axis of the laboratory coordinate frame. 

The derivations of (|4.3|l may be found in |14l Section 5.1] where the authors give 
an argument based on normal modes of phonons. At low enough temperature (i.e., 
the Lamb-Dicke regime), only the lowest order or the CM mode plays a role, and 
this leads to the Hamiltonian 14.3(1 used by Cirac and Zoller in 7 . 

Remark 4.1. It can be shown that the Rabi frequency fin for the internal 
qubit transition coupled by a and a) to the vibration is r]VL/(2\fN), much smaller 
than that, f2, in ((3.1 1|) for a single qubit (rotation) operation. □ 

Lemma 4.1. Let Hj be given by (14.31) for j = 1, 2, . . . , N. Then Hj has a family 
of invariant two-dimensional subspaces 

(4.6) Vj, k = span{\0)j\k + l)cM,\l)j\k)cM}, fc = 0,1,2,..., 
in the Hilbert space 

(4.7) Vj = span{\0) 3 \k) C M,\l)j\k)cM\ fc = 0, 1, 2, . . . }. 

/ oo \-L 

On the orthogonal complement of all the Vj : k 's, i.e., Vj — span [j V^k ) cFj, 
the action of Hj is : Hj \ tp) = for all %j) € Vj . 
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Proof. We have 



H 3 \0) 3 \k + 1) CM 
(4.8) 

and, similarly, 
(4.9) 



2VN 



ll^OlO^-e-^Hfc + 1) C m) + |0),(l|0),e^(at|fc + 1) CM )] 



2VN 



e-^Vk + l |1)#>cm, 



#i|l>#>CM 



hr/Q 



2VN 



+ 1 |0),-|fc + l) C M. 



Thus, Vj,k is an invariant two-dimensional subspace of Hj in Vj. 



□ 



Let the time-evolution operator of Hj (depending on <j)) on Vj be Uj(t,(f>) = 
e~i ff i*. Then on the subspace Vj^, k = 0, 1, 2, . . . , oo, with respect to the ordered 
basis {|0)j|£; + 1)cMj |1)j|^)cm}j Hj admits the following matrix representation 



(4.10) 



HrjQ k + 1 



e- u 




Thus, with Hj restricted to Vj^k, its time-evolution operator is given by 
U hk (t^) = U 3 (tA) 



(4.11) 

Also, note that 
(4.12) 



cos(£fet) — ie ^sm(£kt) 
-ie 1 ^ sin(£fct) cos(£fci) 



h-qVt jk + 1 
~N~' 



Uj(t,<p)\i>) = \i>) for all |V) €Vj, 



because the action of Hj on Vj is annihilation. In physics, these states in Vj are 
said to be off-resonance. 
Next, define 



(4.13) 



2VN 



[H^-e-^a + lO^e^at], j = l,2, 



where |a) is the auxiliary quantum state as indicated in Fig. 14. 2| Then, similarly 
to Lemma 14.11 we have the following. 

Lemma 4.2. Let H* ux be given by (|4"T3|l for j = 1, 2, . . . , JV. Then H? ux has a 
family of invariant two-dimensional subspaces 

(4.14) VJT = 8P*n{\0)j\k + l)cM,\<x)j\k)cMh k = 0,1,2, ... , 
in the Hilbert space 

(4.15) V// ux = s P an{\0)j\k)cM, |a)#>CM \ k = 0,1,2,.. .}. 



On the orthogonal complement V^ aux = 
m ux is 0: m m \4i) = for all S V 3 



span [j V a ^ x I C V 3m , the action of 

k=Q J 

□ 
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The time-evolution operator corresponding to £f aux on V 7 aux is denoted as 
[/ aux (i, </>) = e^^i Now, restrict if aux to the invariant two-dimensional sub- 
space V^ x with ordered basis {|0)j|fc -I- 1)cm, \ot)j |&)cm}; its evolution operator 
has the matrix representation 

t/jXM) = c/; ux (i,^)| yaux = e -* H r*| 

3,k j,k 

_ HrjCt jk + 1 

' £k = — y—- 

Here, we again also have 

(4.17) (7 aux (t,</>)l^) = \i>) for all |^) e F/ ux , 

because ff aux annihilates the subspace Vj aux . Using the CM mode as the bus, we 
can now derive the following 2-bit quantum phase gate. 

Theorem 4.3. Let Uj(t,<f)),Uj tk (t,(t>),U% n *(t,<j>) and U^(t,<j>) be defined as 
above satisfying (|4. 1 1|> <|4. 1T|) for j = 1, 2 and k = 0. Then for 



(4.18) U = J7i(T,0)[/ 2 aux (2T,0)C7i(T,0), T = — Vn, 

rjil 

we have 

(4.19) C/|0) 1 |0) 2 |0) CM = |0)i|0) 2 |0) CM , 

(4.20) C/|l)i|0) 2 |0) C M = |l)i|0) 2 |0) CM , 

(4.21) ^|0) 1 |l) 2 |0)cM = |0)i|l) 2 |0)aM, 

(4.22) U\ 1) x | l)i|0) cm = -|l)i|l)i|0) CM - 

Consequently, by ignoring the last CM-bit \0)cm, we have the phase gate 

(4.23) u = Qtt, cf. until) . 



PROOF. We first verify l|4.19|l : 

c/|o)i|o) 2 |o) CM = f/i(r,o){/ 2 aux (2T,o)[c/i(r,o)|o)i|o)cM]|o) 2 

= U^T, 0)f/ 2 aux (2T,0)[|0) 1 |0) CM ]|0) 2 

(by KWi because |0)i|0) CM e Vj for j = 1) 
= C/ 1 (T,0)[(7 2 aux (2r,0)|0) 2 |0)cM]|0) 1 

- C/i(T,0)[|0> 2 |0) CM ]|0)i 

(by ETfll because |0) 2 |0) CM e V, aux for j = 2) 
= [(7i(T,0)|0)i|0) C m]|0) 2 

- |0)i|0) C m|0) 2 

(by KWi because |0) 2 |0) CM € ^ for j = 1) 
(4.24) = |0) 1 |0) 2 |0)cm- 



(4.16) 



COs(£fci) 



-4p-i<f> 



sin(£ k t) 



-ie 1 ® sin(£ fe i) cos(£ fc i) 



MATHEMATICAL MODELS 



25 



Next, we verify (|4.2()|) : 

t/|i)i|o) 2 |o) CM = [/i(T,o)c/ 2 aux (2r,o)[[/ 1 (r,o)|i)o|o) CM ]|o) 2 

= (7 1 (T,0)C/ 2 aux (2T 1 ,0)H|0)i|l)cM]|0) 2 

(using lUlHl with k = 0, = and £ k T = tt/2) 

= IMT, 0){H)[C/ 2 (2T,0)|0) 2 |1) C m]|0)i} 

= C/ 1 (T,0){[i|0) 2 |l) C M]|0)i} 

(using (|CT3|) and jUJl with ft = 0, <j> = and 2£ fe T = tt) 

= i[CTi(r,0)|0) 1 |l)cM]|0)a 

= ;[H)|i>i|o>cm]|o> 2 

(using @JT) with k = 0, = and £ fc T = tt/2) 

(4.25) = |l)i|0) 2 |0) CM . 

The verification of (|4.21|) can be done in the same way as that for l|4.19|l . 
Finally, we verify (|4.22l) : 

C/|1)i|1) 2 |0) C m - C/i(T,0)C/ 2 aux (2r,0)[C/ 1 (T,0)|l) 1 |0) CM ]|l)2 
= U t (T, 0)C/ 2 aux (2T,0)H|0)i|l) C M]|l> 2 

(using Iginil with j = 1, k = 0, (j) = and £ fc T = tt/2) 
= ^(T, 0)[(-z)^ 2 aux (2T, 0)|1> 2 |1}cm]|0> 1 
= ^i(T,0)[(-i)|l) 2 |l)c M ]|0)i 

(by gJJ because |1) 2 |1)cm G T/, aux for j = 2) 

= (-*) [t^l (r, 0) 1 0> 1 1 1> C!M] 1 1>2 

= H)[(-Oli)i|o)oM]|i) a 

(using gjjl with j = 1, k = 0, <j) = and £ fe T = tt/2) 

(4.26) = -|l)i|l) 2 |0) CM . 

The verifications are complete. □ 

The CNOT gate can now be obtained according to Corollary 12.31 As a conse- 
quence, we have the following. 

Theorem 4.4. The quantum computer made of confined ions in a trap is uni- 
versal. 

PROOF. Use (l2~T5j) and (|4"2"3|) to deduce Corollary IO □ 

A very similar gate has been implemented experimentally using a cold, confined 
Bc + ion |19| . To conclude this section, we mention two more schemes using ion 
traps. 

(a) The S0rensen-M0hner scheme 26 . 

Experimental inconveniences in implementation of the Cirac-Zoller scheme on con- 
fined ions include the requirement that the ions be initialized in the ground vi- 
brational state |0)cm, which requires sophisticated laser cooling techniques, and 
the experimental observation that the CM vibration of a linear string of ions in 
a trap has limited coherence. The S0rensen-M0lmer scheme applies to a string 
of ions, but the ions can have thermal motion, as long as they remain within the 
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Lamb-Dicke regime. The ions never populate a vibrational mode, so decoherence 
effects are reduced. The idea is to simultaneously excite an ion pair using two laser 
beams detuned by frequency S from the upper and lower vibrational sidebands 
u>o + wqi — S and loq — wcm + <5, so that the energy (or frequency) sum is 2u>q, 
corresponding to excitation of both ions by two interfering transition paths from 
|0, 0) \ri) to |1, 1) \n), without changing the vibrational state. The lasers are detuned 
from the red and blue sideband transitions of single ion excitation, so the vibration 
is only virtually (negligibly) excited, and the vibration levels are never actually pop- 
ulated. One laser beam is applied to each ion. For any two ions, each illuminated 
by both of the laser beams, an additional two interference paths exist, and the effec- 
tive Rabi frequency for the two-ion transition becomes f2sM = -(i)fi) 2 /( w CM — S), 
where f2 represents the single ion Rabi frequency at the same laser power. The 
two-ion interaction Hamiltonian can be written (M2sm/2)(|0, 0) (1, 1| + |1, 1)(0,0| + 
|1,0)(0,1| + |0,1)(1,0|). The Rabi frequency can be approximated as Qsm = rfVL 2 /S, 
since wcm — 5 ~ 5. Full entanglement is achieved at a particular interaction time 
T = 7r/2f2sM, when the wavefunction becomes a coherent superposition of max- 
imally entangled |0,0) and |1, 1) states \ip) = 2" 1 / 2 (|0,0) -i|l,l)). Once entan- 
glement is achieved, the CNOT gate can be implemented with additional single 
particle rotations, carried out analogously to those in the Cirac-Zoller scheme, 
with individual laser beams exciting single ions. Simultaneous entanglement of 
any even number of ions can be carried out similarly, although the CNOT gate 
applied to any pair of ions is sufficient for computation. The technique has been 
experimentally demonstrated with four ions in a small linear ion trap 22 . 



(b) The Jonathan-Plenio scheme |15| . 

When an atomic transition is driven by a laser beam at resonance, an AC Stark shift 
(or light shift) results, changing the energies of the coupled levels due to the coupling 
to the laser field. The laser drive on an ion in a linear string in a trap produces 
Rabi transitions between the coupled levels, but if the Rabi frequency matches one 
of the CM vibration modes v v of the ion string, an exchange of excitation energy 
between the internal and vibrational variables also occurs. This exchange can be 
used to derive a CNOT gate between any two ions in the linear chain, based on 
an interference analogous to that of the S0rensen-M0lmer scheme discussed above. 
This is expected to result in relatively fast quantum gates, but requires cold ions. 
When the light shift is used to drive virtual two-phonon transitions simultaneously, 
using more than one vibrational mode, it is predicted that gate speeds will be 
somewhat reduced, but the ions remain insensitive to heating, possibly offering an 
advantage over the S0rensen-M0lmer scheme. Two ion states of the form |1, 0)|n p ) 
and |0, 1)1%,) are coupled. The effective time-independent Hamiltonian is written 
— 0)(0, 1| + |0, 1) (1 , 0|] plus a sum over states p which depend on the ion 
motion, but which can be cancelled by reversing the sign of a parameter half way 
through the operation. The states in the Hamiltonian are assumed now to be 
"dressed", i.e., formed from the original states coupled and shifted by interaction 
with the external field. The frequency is u — (tt 2 / tyYip^rjip^pVp / {Q 2 — v 2 )- 
Maximally entangled states can be produced at the time T — |7r/4w|. This proposed 
technique has not yet been experimentally demonstrated. 
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5. Quantum Dots 



Quantum dots are fabricated from semiconductor materials, metals, or small 
molecules. They work by confining electric charge quanta (i.e., spins) in three 
dimensional boxes with electrostatic potentials. The spin of a charge quantum in a 
single quantum dot can be manipulated, i.e., single qubit operations, by applying 
pulsed local electromagnetic fields, through a scanning-probe tip, for example. Two- 
qubit operations can be achieved by spectroscopic manipulation or by a purely 
electrical gating of the tunneling barrier between neighboring quantum dots. 

The spin-1/2 operator of an electron is given by 



(5.1) 



where a x , a y and a z 



are the usual Pauli matrices: 



(5.2) 
and 

(5.3) 








r 




"0 -i 




'1 " 




1 







i 


, <r z = 


-1 



"l" 




"o" 




"0" 





, Gy — 


1 


, e 2 = 















1 



are the unit vectors in the directions of x,y and z. Let Si and Sj denote, re- 
spectively, the spin of the electric charge quanta at, respectively, the i-th and j-th 
location of the quantum dots. Then the usual physics of the Hubbard model 
gives the Hamiltonian of the system of n quantum dots as 



(5.4) 



H = J2^B9 j {t)B j (t)-S j 



^2 J jk{t)Sj 

l<j<k<n 



where the first summation denotes the sum of energy due to the application of a 
magnetic field Bj to the electron spin at dot j, while the second denotes the inter- 
action Hamiltonian through the tunneling effect of a gate voltage applied between 
the dots, and 



Jjk (t) 



(5.5) 

Note that for 

(5.6) 

and 



is the Bohr magneton; 

is the effective g-i actor; 

is the applied magnetic field; 

the time-dependent exchange constant [181 see [10] in the Refer- 
ences therein], with Jjk{t) = 4^ fe (t)/u, which is produced by the 
turning on and off of the tunneling matrix element tij (t) between 
quantum dots i and j, with u being the charging energy of a single 
dot. 



■ a y 3) e y + cr«e 2 



i = 1,2,. 



(5.7) Bj{t) = b^(t)e x + b^(t)e y + b^(t)e z , j = 1, 2, 
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the dot products in Ij5.4(l are defined by 

(5.8) B 3 (t) ■ Si = bP(t)*P + bV\t)<r<p + b¥\t)<rP. 

In Fig. 15. II we include a quantum dot array design given in Burkard, Engel and 
Loss |6]. 

Again, by the universal quantum computing theorems in §2, we need only 
discuss the case of a system with two quantum dots. The spin-up f and spin-down 
I of the electric charge quanta in each dot are denoted as |0) and |1), respectively. 
Then the underlying Hilbert space is 

(5.9) W = B P «m{|00),|01),|10),|ll)}. 
The wave function \ip(t)) of the system has four components: 

(5.10) m)) = M*)|00) + ^(t)|01) + V 3 (*)|10) + ^(t)|ll>, 
satisfying the Schrodinger equation 

(511) f ih^ t \m)=H(t)m)), t>o, 

\ \ip(o)) = \ipo) e H, 

where 

(5.12) H{t) = ^[n 1 (t)-a+Sl 2 {t)-T + Lj(t)a-T], 
following from (|5.4|l by rewriting the notation 

(5.13) S 1= a, S 2 =t; i i B g j {t)B j {t) = hl j {t), j = l,2; J 12 (t) = |w(t). 

The Qi(t),Sl2(t) and u>(t) in (|5.12|l are the control pulses. 

Let us now give a brief heuristic physics discussion for the derivation of the 
Hamiltonian in (|B .4|> and Q5.12JI as a remark below. 



Remark 5.1. The suggestion to use quantum dots as qubits for the purposes of 
quantum information processing in general, and quantum computing in particular, 
is based on the following observations. To set the stage, we note that a quantum 
dot is formed by a single electron that is tightly bound to a "center" - some dis- 
tinguished small area in a solid. Typically, the binding is strongly confining in one 
spatial direction (the z direction) while the electron has some freedom of motion 
in the plane perpendicular to it (the xy plane), with its dynamics governed by 
the forces that bind it to its center. In addition, there is the possibility of ap- 
plying electric and magnetic fields of varying strength and direction, whereby the 
experimenter can exert some external control over the evolution of the electron's 
state. 

In particular, a magnetic field B(t) will couple to the magnetic moment of the 
electron, which gives rise to a term in the Hamilton operator that is proportional 
to B ■ cr, where a is the Pauli spin vector operator of the electron, so that 

(5.14) H = -Ml(t)-a 
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Figure 5.1. Quantum dot array, controlled by electrical gating. 
The electrodes (dark gray) define quantum dots (circles) by con- 
fining electrons. The spin 1/2 ground state (arrow) of the dot 
represents the qubit. These electrons can be moved by electrical 
gating into the magnetized or high-g layer, producing locally dif- 
ferent Zeeman splittings. Alternatively, magnetic field gradients 
can be applied, as e.g. produced by a current wire (indicated on 
the left of the dot-array). Then, since every dot-spin is subjected 
to a different Zeeman splitting, the spins can be addressed indi- 
vidually, e.g. through ESR (electron spin resonance) pulses of an 
additional in-plane magnetic ac field with the corresponding Lar- 
mor frequency lul = qfJ*BB±/h. (In the figure, B± denotes the 
component of the magnetic field perpendicular to the array plane, 
while B\ | denotes that parallel to the plane.) Such mechanisms can 
be used for single-spin rotations and the initialization step. The 
exchange coupling between the quantum dots can be controlled by 
lowering the tunnel barrier between the dots. In this figure, the two 
rightmost dots are drawn schematically as tunnel- coupled. Such 
an exchange mechanism can be used for the CNOT-gate operation 
involving two nearest neighbor qubits. The CNOT operator be- 
tween distant qubits is achieved by swapping (via exchange) the 
qubits first to a nearest neighbor position. The read-out of the 
spin state can be achieved via spin-dependent tunneling and SET 
devices |18) . or via a transport current passing the dot |18| . Note 
that all spin operations, single and two spin operations, and spin 
read-out, are controlled electrically via the charge of the electron 
and not via the magnetic moment of the spin. Thus, no control of 
local magnetic fields is required, and the spin is only used for stor- 
ing the information. This spin-to-charge conversion is based on 
the Pauli principle and Coulomb interaction and allows for very 
fast switching times (typically picoseconds). A further advantage 
of this scheme is its scalability into an array of arbitrary size. (The 
figure is excerpted from Burkard, Engel and Loss [§], courtesy of 
Springer- Verlag Berlin Heidelberg, while the caption is excerpted 
from Schliemann and Loss L 23 .) 
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is the structure of the effective single-qubit Hamilton operator, with (l(t) oc B(t) 
being externally controlled by the experimenter. By suitably choosing B(t) various 
single-qubit gates can be realized. 

To suppress the unwanted effect of the Lorentz force (which couples to the 
electron's charge, results in a time-dependent energy, and thus produces an accu- 
mulated phase change), it is expedient to have the magnetic field vector in the xy 
plane, so that B ■ e z — is imposed. Nevertheless, general single qubit gates can be 
realized by a succession of rotations of a about axes in the xy plane. This is just 
an application of Euler's classical result that general rotations can be composed of 
successive rotations about two orthogonal axes only (standard choice is the z axis 
and the x axis but, of course, the x axis and the y axis serve this purpose just as 
well). 

More complicated is the implementation of two-qubit gates. Here one needs two 
quantum dots in close proximity, by which is meant that they have a considerable 
interaction, predominantly originating in the Coulomb repulsion of the charges of 
the two electrons. It is important that external electric fields can be used to modify 
the effective interaction strength in a wide range, and in particular one can "turn 
it off" by separating the two quantum dots (or rather shielding them from each 
other). 

Consider, thus, two quantum dots with some interaction potential in addition 
to the potentials that bind them to their respective centers. This situation is 
reminiscent of the hydrogen dimer, the H2 molecule, except that the confinement to 
a plane and the form of the binding potential, and also of the effective interaction, 
are quite different. It is, indeed, not so simple to model the various potentials 
reasonably well, and one must be practiced in the art of solid-state theory to do it 
well. We shall, therefore, refer the interested reader to the specialized literature, 
perhaps starting with p>] and the references therein. 

The general picture, however, can be grasped without getting involved with 
such details. First, we note that at sufficiently low temperatures, only the ground 
state and the first excited state of the interacting two-dot system will be dynam- 
ically relevant, at least as long as we carefully avoid exciting other states. For 
symmetry reasons, the ground state is a spin-singlet (total spin angular momentum 
of Oh) and has a symmetric spatial wave function, whereas the first excited state is a 
spin-triplet (total spin angular momentum of lh) and has an antisymmetric spatial 
wave function. The excited state is long-lived because triplet-to-singlet transitions 
tend to have very small matrix elements. 

The total spin angular momentum vector operator is S = \h(a + r) where we 
denote the Pauli vector operators of the two electrons by a and r, respectively. The 
eigenvalues of S 2 are Oh 2 in the spin singlet and 2h 2 in the spin triplet. In view of 
a 2 = 31 and r 2 =31, this says that a ■ t has eigenvalue —3 in the singlet and +1 
in the triplet. As a consequence, 



Effectively, then, they project on the ground state and the excited state, respec- 
tively. Denoting by E the ground state energy, and by E\ = E + 2hw that of the 



(5.15) 
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excited state, we have 

H = ^(1 - a ■ t)E + -(31 + a ■ t)E! 

(5.16) =^(E + 3E 1 )1 + \(E 1 -E )a-T 
or, after dropping the irrelevant additive constant, 

(5.17) H = -fvjo{t)o ■ t 

for the effective Hamilton operator of the coupling between the two quantum dots. 
We have indicated the externally controlled time dependence of u>, the "turning 
on" and "turning off" of the interaction at the experimenter's discretion, by t. 

In summary, then, we have two terms of the form 1)5. 14[) . one for each quantum 
dot, and the interaction term l|5.17|l in the effective Hamilton operator 

(5.18) H = -h[Qi(t) ■ a + n 2 (i) • r + u(t)a ■ r] . 

Within reasonable ranges, the experimenter is capable of realizing any ili(i), ^2(^)1 
and wit), where it is fully sufficient to have two of them vanishing at any instant. 
□ 



Lemma 5.1. Let a and r be defined as in (|5.13[) . Then we have 
(5-19) (i) f— =1 and (a -t) 2 = '31-2<t -T. 

(5.20) (ii) U sw = — (1 +<r • r) = U^ w = U^, where U sw is the swapping gate 

(5.21) U sw \jk) = \kj) for j, k = 0, 1; cf. W^W - 
(iii) 

(5.22) UlaU sw = r, U^tU sw = r 



(5.23) (iv) C/ S 2 W = 1. 

Proof. Let us recall first how the various <jj and Tj transform the basis vectors, 
for j = x,y,z: 





^|o-) = |i-; 


>, 




<7x|l-) = |0-), 




<Ty\0') =*|1 


•>, 




o- y \l-) = -i|0- 


(5.24) 


^|o-) = |o-; 


>, 




^|l-) = -|l-) 


and likewise for 


t x , T y , and 


t z acting on | • 


0) and 


|-1). Thus 




a -t|00) = 


= (1 2 + * 2 )|11> 


+ |00) = 


= 100} , 




ct-t|01) = 


= (l 2 -z 2 )|10) 


-|01> = 


= 2|10)-|01), 




ct-t|10) = 


= (l 2 -z 2 )|01) 


+ |10) = 


= 2|01> - |10> , 


(5.25) 


a-T\ll) -- 


= (l 2 + z 2 )|00) 


+ 111) = 


- Ill) , 
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and we see that the 4x4 matrix for a ■ r is given by 

"1 0' 



(5.26) 



a t = 



0-12 
2-1 




and the first equation in Ij5.19|) follows. 

We read off that a ■ r has a 3- fold eigenvalue +1 and a single eigenvalue —3. 
The respective projectors to the subspaces characterized by these eigenvalues are 
easily verified to be 



(5.27) 



-(31+a-r) and P 2 = -(1 - a -t) , 



where P 2 



4' 

"j for j — 1, 2, and PjPfc = for j ^ k. Therefore, for any sufficiently 
well-behaved function, including polynomials and analytic functions, of a -t, 

f(a ■ t) = /(l)i(31 + a ■ t) + /(-3)i(l -a-r) 



(5.28) 



= \ [3/(1) + /(-3)]l + i[/(l)-/(-3)] <t t. 



according to the spectral theorem. As a first application, consider f{x) 
find 



= x 2 and 



(5.29) 



(a ■ t) 2 = 31 - 2a • r , 



which is the second equation in H5.19|l . 
Further, in view of the matrix for U sw , 



(5.30) 



a • r ) = 



10 
10 
10 
1 



the swapping gate really swaps: U sw \jk) — \kj) (no effect on |00) and |11), whereas 
1 01} and 1 10) are interchanged). To verify 

(5.31) U sw cr = tU sw 



it should suffice to inspect the pair 
(5.32) u x = |00)<10| + |01)(11| 

and 



|10)(00| + |11)(01| 



(5.33) 



= |00)(01| + |10)(11| + |01)(00| + |n)(io| 



for example. One can also have a purely algebraic, but more lengthy, argument 
that exploits nothing but the basic identities, such as a 2 = 1 and a\a2 — icr^, etc. 
Thus, the rest easily follows. □ 

Theorem 5.2. Denote by U{t) the time evolution operator for the quantum 
system l|5.11|l and l|5.12|) for time duration t £ [0, T]. Choose fti(t) = Cl2{t) = in 
(|5.12|l and let to(t) therein satisfies 



(5.34) 



u(t)dt 



IT 

2' 
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Then we have U(T) = — e 7!l / i U sw , i.e., U(T) is the swapping gate (with a nonessen- 
tial phase factor —e nt / i .) 

Proof. By assumptions, we have now 

(5.35) H(t)=w(t)a-T/2. 
Since u(t) is scalar- valued, we have the commutativity 

(5.36) ff(ti)fT(ia) = H{t 2 )H(tx), for any t 1: h 6 [0,T]. 
Therefore 

U(T) = e - i /o TiT (*) dt /fi = e [~lf/o T «(*)*]»-T 



1 " T 



uj(t)dt 



(5.37) = cos(0ct • t) — isin(<^<7 ■ t), 



where e l ^* 7 ' T ; cos(</xt • t) and singer • r) are 4x4 matrices. We now utilize (|5.28|) 
to calculate the exponential matrix U(T): 

(5.38) U(T) = e-^ a - T = e'** ■ |(31 + a • r) + e~ 3 ^ ■ ^(1 - a ■ r). 

Thus, with a little manipulation, (|5.38(l becomes 

1 +<7 - T 



t/(T) = e 



cos(20)l — isin(20)- 



2 

(5.39) = e^[cos(2<£)l - i sin(20)[/ sw ], by f5~2U|) . 

Choosing = 7r/4, we obtain the desired conclusion. □ 

1/2 

Corollary 5.3. The square roots of the swapping gate, U S w , are 

±TTi/A 

(5.40) UV? = — ^-(lT iU sw ). 

PROOF. From l|5.39() . we first obtain 

(5.41) U sw =ie"T[/(T). 
Then use = ±7r/8 in (|5.39|l to obtain 

and the desired conclusion. (Note that these two square roots of U sw reflect the 
choices of \/l — 1 and 

(5.43) the square root of —1 = ±i 

for the square roots of the eigenvalues of U sw .) □ 

Theorem 15.21 gives us the choice of control pulses ui(t) (with fii(t) and 02 (t) 
being set to 0) which yields the swapping gate J7 SW , which does not cause entan- 
glement. However, C/ S w 2 in Corollary 15 . 31 causes entanglement. See Corollary 15.51 
shortly. 



(5.42) Uli 2 = (ie-T)V2 e ±«/s 



34 



G. CHEN, D. A. CHURCH, B.-G. ENGLERT, AND M. S. ZUBAIRY 



Theorem 5.4. Letcf>,9 G [0, 2ir] be given. Denote e(0) = cos <f>e x + sin <fie y + 0e z 
for the given (f>. Let Ui.OiCO be the time evolution operator corresponding to the 
quantum system (|5.11|l and <|5.12[1 for t £ [0, T] where the pulses are chosen such 
that 



(5.44) Qx(t) =fii(i)e(0), fi 2 (*)=0, w(t) = 0, 
with Qi(t) satisfying 



te[o,T], 



(5.45) 



S7i (t)c?t = 20, /or £/ie given 



Then the action of U\^ t (t) on the first qubit satisfies 
(5.46) Ui^i(t) = E/g.0, i/ie unitary rotation gate \2.17]) . 

Proof. We have 



cos f 
-ie^ sinf 



COS ( 



cos 6>1 



ie -i * sinf 



ie 10 sin* 



(5.47) 



cos 91 — i sin 9 cos 0(7^ — i sin sin <pa y 
cos 91 — i sin ^(cos ^cr^ + sin </>o- y ) 
cos 91 — i sin 6*e(</>) • a 

-i9e(<p)-a 



noting that in the above, we have utilized the fact that the 2x2 matrix 



(5.48) 



e(4>) a 







COS ( 



I sin ( 



cos < 



I sin < 







satisfies (e(0) • a) 2n = 1 for n = 0, 1, 2, . . . . 

With the choices of the pulses as given in (|5.44|) . we see that the second qubit 
remains steady in the time-evolution of the system. The Hamiltonian, now, is 



(5.49) 



#i(i) = -fii(i)ei(0)-ff 



and acts only on the first qubit (where the subscript 1 oie\{4>) denotes that this is 
the vector e{<f) for the first bit). Because f2i(t) is scalar- valued, we have 



(5.50) 
Thus 



(5.51) 



#l(tl).ffl(t2 



Hi(t 2 )ffi(*i) for any h,t 2 G [0,T]. 



e-i/o T °i(*) e i(^)-' T dt 
e [-i/o rfi i(*) d *] e i(</')- 0. 



e 



(by 



using (|5.47|l . The proof is complete. 

We may define L/2,n 2 m a similar way as in Theorem 15.41 



□ 
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Corollary 5.5. The quantum phase gate (|2.18|) is given by 
(5.52) Qi, = (-i)U hn?) U 2 ,a 2 U^ 2 U hn wU^ 2 , 

where 

J Sl^\t)dt=-Tre lz , 
(5.58) { J il ( i } (t)dt = ne lz /2, 

J n 2 (t) dt = -ne 2z /2, 

and e\ z ,e2z denote the e z vector of, respectively, the first and the second qubit. 

Remark 5.2. In order to realize this succession of gates, only one of the Sl(t) 

in l|5.58|l is nonzero at any given instant t, with the duration when Sl^ (t) 

earlier than that when (I2W 7^ 0, and that when Cl\ (t) ^ even later. Earliest is 

1/2 

the period when ui(t) ^ for the first U B w , and another period when u>(t) ^ is 
intermediate between those when Oj^i) ^ and fi2(t) 7^ 0. 

Proof. Define 

(5.54) Uxor = e^'e^U^e^U^ 2 , 

1/2 — — i 

with U S w — e Jt (l + iU svr ) chosen from l|5.4U|l . Then it is straightforward to check 
that 

EW|00) = |00)(i), C/xor|01) - loixo, 

(5.55) CW|10> = |10)(i), Dxor|H> = |H)(-0. 
so that 

Uxor = f(|00><00| + |01)(01| + |10)(10| - |11>(11|) 

(5.56) = iQ„. 

□ 

Remark 5.3. In |18l eq. (2)], Loss and DiVincenzo write their Uxor gate as 

(5.57) Uxor = e^'e^'U^^'uU 2 . 

Thus there is a difference in notation between the above and (|5.54|) : 

(5.58) Si = a z /2, S| = t z /2. 

□ 

Theorem 5.6. The quantum computer made of quantum dots is universal. 

PROOF. This follows easily from l|2.19|) . Theorem l5.4l Corollarv l5.5l and finally 
Corollary IO □ 
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Throughout this paper we have been employing the Schrodinger point of view 
and could therefore afford to identify the ket vectors with particular numerical 
column vectors, the bra vectors with the corresponding row vectors, and we did 
not distinguish between operators and matrices. Now, let us give a brief, heuristic 
exposition from the point of view of the Heisenberg picture, in which it is important 
to discriminate between kets and their columns, bras and their rows, operators and 
their matrices. Dirac's classic text and the recent textbook by Schwinger |24| 
treat quantum mechanics consistently in the Heisenberg picture. 

For a start we return to l|5.9l) and (|5.1U|) where each \jk) is a joint eigenket 
of o~ z and t z with eigenvalues 1 — 2j and 1 — 2k, respectively for j, k = 0,1. A 
simultaneous measurement of o~ z and r 2 yields the actual one of the four pairs of 
values and tells us which of the four possibilities is the case. 

Now, a measurement of o~ z , say, at time t\ is not the same as a measurement 
of o z at time £2, one really has to perform two different experiments (perhaps 
using the same lab equipment twice). Therefore, according to the reasoning of 
Heisenberg and Dirac, just to mention the main quantum physicists, we should also 
have two mathematical symbols for these two measurements: o z (t\) and u z (t-2). 
As a consequence, then, also the eigenkets must refer to the time of measurement: 
\jk,ti) and \jk, ti). Accordingly, the eigenvector equations read 

a-Jt)\jk,t) = \jk.t)(l - 2j) ) 

(5.59) W ' ' ' A ' UoralH, for j.k = 0,l. 
r z (t)\jk,t) = \jk,t)(l-2k) J J ' 

On the other hand, the state of the system, which is typically specified by a state- 
ment of the kind "at time to both spins were up in the z direction" (formally: state 
ket I VP) = |00, to)), does not depend on time by its nature. So, rather than JOJ5J, 
we write 

(5.60) I*) =^ 1 (t)|00,t)+^(t)|01,t)+^ 3 (t)|10,t)+^(*)|ll,*) ) 

where the dynamical time dependences of the kets |00, t), . . . , |11, t) is compensated 
for by the parametric time dependence of the probability amplitudes i>i(t), . . . , 
"4>4,{t), so that the state ket \^>) is constant in time, as it should be. 
The probability amplitudes are thus given by 

(5.61) 

^i(i) = <oo,t|*), ^(i) = (0M|*), = <i<M|*> , ik(t) = <iM|*), 

and they have the same meaning as before: |-02(^)| 2 , for example, is the probability 
that the first spin is found up and the second down at time t. And, of course, they 
obey the same Schrodinger equation as before, namely (|B . 1 1|) where \^{t)) is just 
the column [ipi(t) ip2(t) ip3(t) ipi(t)] and H(t) the 4x4 matrix that is meant in 

With the Hamilton operator (not matrix!) 

(5.62) H(a(t),r{t),t) = ^[n 1 (t)-a(t)+a 2 (t)-T(t)+u(t)a(t)-T(t)}, 

which has a dynamical time dependence in a(t) and r(t) and a parametric time 
dependence in ft^(t), j — 1,2, and Lo(t), the four rows of this Schrodinger equation 
are just 

(5.63) ih^(jk,t\V) = (jk,t\H(v(t),T(t),t)\V) for jfc = 00, 01, 10, 11 
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and the elements of the 4x4 matrix of (|5.12(l are 
(5.64) 

(jk,t\H(a(t),T(t),t)\lm,t) with jk = 00,01, 10, 11 and Im = 00,01, 10, 11. 

The Pauli matrices in (|5.2|l appear when we express the components of <j(t) 
explicitly in terms of the eigenkets and eigenbras of cr z (t) and r z (t), 



fc=0,l 



E (|0fe,t)(lfc,i| + \lk,t){0k,t\) 

[(0k,t\,(lk,t\] 



(5.65) 



(7y(t) 

o z (i) 



E 

fc=0,l 

E 

fc=0,l 

E 

fc=0,l 



0fc,i) 
1/M) 



1 

1 



0k,t) 
lk,t) 

0k,t) 
lk,t) 



[<0*,t|,<lfe,*|] , 

[<0M|,<iMI] 



'\jo,t)~ 




'0 1" 




'0 -i 




'1 " 




\jl,t)_ 


( 


1 


e x + 


i 


e y + 


-1 





[0'o,*l,0'i.*l] 



and likewise for r(i), compactly written as 
(5.66) 

r(t) = E 

3=0,1 

Note how numerical 2x2 matrices are sandwiched by 2-component columns of kets 
and 2-component rows of bras to form linear combinations of products of the dyadic 
ket-times-bra type. 

The Schrodinger equation l|5.63|l is the equation of motion of the bras (jk,t\ 
and those of the kets \jk, t) is the adjoint system. The equation of motion for the 
operator a{t) follows from them, 



(5.67) 
and likewise for r(t), 
(5.68) 



ih-a(t) = a(t)H t - H t a{t) - [a(t),H t ] , 



ih-T(t)=[T(t),H t ], 



where H t abbreviates H(a(t),T(t),t) , the Hamilton operator of (|5.62() . These are 
particular examples of the more general Heisenberg equation of motion, which states 
that any operator - here: any function Ot of a(t), r(t), and t itself - obeys 



(5.69) 



where ^ is the total time derivative whereas differentiates only the parametric 
time derivative, so that the commutator term accounts for the dynamical change 
in time of Ot- Note that (|5.69|l has exactly the same structure as the Hamilton 
equation of motion in classical mechanics. 

A particular O t is the unitary evolution operator Ut ,t = U(to;cr(t),T(t), t) that 
transforms a bra at time to into the corresponding bra at time t (earlier or later) , 



(5.70) 



to,t i 
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where the ellipsis stands for any corresponding set of quantum numbers (such as 
00, . . . 11). Ut ,t has a parametric dependence on both to and t, in addition to 
the dynamical t dependence that originates in a(t) and r(t). The group property 
Ut„.t 1 Ut lt t 2 = Ut Qt t 2 follows from the definition and so does the unitarity, and thus 

(5-71) U t0 ^ = U to ,t~ l = U t j 

is implied immediately. Accordingly, we have 

(5.72) a(t) = U t .t a{t )Ut ,t , r(f) = U t ,t T{t )U to ,t , 

and therefore U(to;a(t),T(t),t^ — U(t ; cr(t ), r(i ), t) , that is to say: it doesn't 
matter if we regard Ut t t as a function of the spin operators at the initial or the 
final time. 

Upon combining (|5.70|) with the Schrodinger equation (|5.63|) . we find first 

(5.73) i^U to ,t = U t0tt H t , 

at 

and then with Heisenberg's general equation of motion l|5.69ll . 
(5-74) iH §i U *0't = HtUt ^ ■ 

To begin with the implicit a and r are all at time t, but since the derivative in (|5.74(l 
refers to the parametric dependence on t, not to the dynamical one, we can just as 
well take all a and r at time to, so that the t dependence is then only parametric. 

With the understanding, then, that we are only dealing with parametric t 
dependences, we can combine (15.7411 with the initial condition Ut 0t t = 1 into the 
integral equation 

(5.75) u t0 , t = l-- / dhH tl U to , tl , 

J to 

which can be iterated to produce a formal series expansion, a variant of both the 
Born series of scattering theory and the Dyson series of field theory, of Ut ,t m 
powers of the (time-integrated) Hamilton operator. It begins with the terms 

(5.76) U to ,t = 1 - ~ / dhH tl -Lfdtxf 1 dt 2 H tl H t2 + ■■■ . 

h J t V J to J t Q 

For the Hamilton operator l|5.62|l . the Heisenberg equations of motion (|5.67() 
and H5.68J) read 

[ ^-a(t) = n 1 (t)xa(t)-uj(t)a(t)XT(t), 

(5.77) <^ <f 

I — T(t) = fl 2 (t) X T(t) + L0{t)(j{t) x r(t) . 

Example 5.1. Let us give the solution of H5.74JI for the case when Cli(t) = 0, 
fi 2 (i)=0: 

(5.78) | ( | 
<*(t) = 2 [ CT (*o) + T (^o)] + ^ [ CT (*o) - r(t )] cos tp - -a(t Q ) x r(t ) sin tp , 

T~(t) = ^ [o-(to) + r(t )] - i [a(t ) - r(t )] cos<^+ -a(t ) x r(i ) sin <p, 
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where 

(5.79) ip = 2 [ dt'u(t') 



One can verify by inspection that (|5.78|) solves (|5.77l) . □ 
Summarizing the above discussions, we have established the following. 
Theorem 5.7. Let (|5.11fl and (|5.12[) be satisfied. Then we have 



(5.80) 




— \a,H](t) = Slx(t) x ait) - u(t)cr(t) x r(t), 
in 

^r[T,H](i) = il 2 (t) x r(t) + w(t)a(t) XT(t) 
in 



an 



d 



(5 - 81) < r( t j = c/ t ; t0 T(to)c/ to ,; ' /orany <o ' t -°- 



6. Conclusion 



□ 



New or improved versions of the elementary devices cited in this paper are 
constantly emerging. For example, in the design of cavity QED, instead of using a 
particle beam to shoot atoms through an optical cavity as illustrated in Fig. 3.2, 
there have been proposals to enclose ion traps inside a cavity. Thus, such proposals 
incorporate both features of ion traps and cavity QED as studied in this paper. 
Even though it is not possible for us to describe all major types of contemporary 
quantum computing devices or proposals, we hope that our mathematical analysis 
and derivations herein have provided reasonably satisfactory rigor and basic prin- 
ciples useful in helping future interdisciplinary research in quantum computation. 

From the control-theoretic point of view, the best theoretical model for laser 
and quantum control seems to be the Schrodinger equation offered by l|5.11[) with 
the Hamiltonian given by 1(5. 4j) . wherein the Bj(i)'s and Jjk(tys are the bilinear 
controllers. Equivalently, the Heisenberg's equations of motion as presented in 
Theorem 15.71 ( which deal with only two quantum dots) form a bi-trilinear control 
system with f2i(t) and £^2 (^) therein as bilinear controllers, while to(t) is a trilinear 
controller. Such control systems (as presented in Theorem 15.71 e.g.) look nearly 
identical to the classical rigid body rotation dynamics and, thus, a fair portion 
of the existing control results appear to be readily applicable. In our discussion 
in Section we have derived certain formulas for the shaping of laser pulses (see 
(I5.34|> . (|5.44l) . I5.45fl . I|5.53[) . etc.) which can achieve the desired elementary 1-bit 
and 2-bit elementary quantum gates. However, we wish to emphasize that such laser 
pulse shaping is underutilized as far as the control effects are concerned. To obtain 
more optimal control effects, feedback strategies need to be evaluated. Quantum 
feedback control is currently an important topic in quantum science and technology. 

For the models of 2-level atoms, cavity QED and ion traps, our discussions here 
have shown that the admissible class of controllers is quite limited, the reason being 
that the control or excitation laser must operate at specified frequencies in order 
to stay within the two or three levels, and to avoid the excitation of unmodeled 
dynamics. We can exercise control mainly through the duration of the activation 
of laser pulses, such as (13.20(1 has shown. 
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There are many interesting mathematical problems in the study of laser-driven 
control systems. Quantum computing/computer contains an abundant source of 
such problems for mathematicians and control theorists to explore. 
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